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Data collected from a worldwide 1992 experiment have been processed at JPL to 
determine precise orbits for the satellites of the Global Positioning System (GPS). A 
filtering technique has been tested to improve modeling of solar-radiation pressure 
force parameters for GPS satellites. The new approach improves orbit quality for 
eclipsing satellites by a factor of two, with typical results in the 25- to 50-cm range . 
The resultant GPS-based estimates for geocentric coordinates of the tracking sites, 
which include the three DSN sites, are accurate to 2 to 8 cm, roughly equivalent to 
3 to 10 nrad of angular measure. 


I. Introduction 

The satellites of the Global Positioning System (GPS) 
are maintained by the U.S. Department of Defense for 
navigational purposes. These satellites are distributed 
in six evenly spaced orbit planes, at an orbit altitude 
of 20,000 km and with an orbit period of approximately 
12 hr. GPS measurements collected from globally dis- 
tributed ground receivers are also being used to estimate 
Earth rotation, polar motion, and station coordinates. Es- 
timation of the geocenter, or Earth’s center of mass, has 
geophysical and scientific implications. GPS estimates of 
the geocenter can be used for precise reference- frame cal- 
ibration and alignment and to eliminate a bias-type er- 
ror in DSN site coordinates. Over time intervals of weeks 
to months, the GPS measurements can be used to pre- 
cisely monitor variations in tracking site coordinates due to 
crustal motion and continental drift. The data used in this 
analysis are taken from the International Global Position- 
ing System Geodynamics Service 1992 (IGS , 92) campaign. 


IGS’92 consisted of approximately 30 globally distributed 
Rogue receiver sites tracking the 18 GPS satellites active 
during this time. The focus of this analysis was to assess 
the effects of mismodeling satellite force parameters due 
to solar radiation on the estimation of other parameters, 
such as GPS orbits and geocentric DSN coordinates. 

II. Estimation Strategy 

A unique strength of GPS measurements is that the 
satellites are sensitive to the geocenter, yet relatively in- 
sensitive to errors in the gravity field because of their high 
orbit altitude and the relatively short data arcs (30 hr) 
needed for the solution. Mismodeling of satellite force pa- 
rameters, however, can have a significant effect on satel- 
lite orbits, especially in orbit prediction [1]. The goal of 
this study was to improve satellite modeling in order to 
achieve centimeter-level accuracy for global geocentric co- 
ordinates. To achieve this level of accuracy, it is essential 
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to include corrections for numerous potential errors, in- 
cluding Earth’s rotation and orientation, atmospheric dis- 
tortion of the radio signals from the satellites, gravitational 
and nongravitational forces acting on the satellites, and 
various geophysical effects. 


technique has enabled the achievement of few-centi meter 
geocentric coordinate accuracy [6]. 

III. Results and Discussion 


The data used in this analysis were taken from GPS 
week 660, which ran from August 30 through Septem- 
ber 5, 1992. In general, the data contain carrier phase 
and pseudo-range measurements from 18 available GPS 
satellites tracked by approximately 30 globally distributed 
JPL Rogue receivers. These tracking sites are shown in 
Fig. 1 and listed in Table 1. The data were processed 
using the GPS Inferred Positioning SYstem/Orbit Analy- 
sis and Software Integration System, (GIPSY/OASIS II) 
software [2,3]. 1 All nonfiducial station locations were esti- 
mated, as well as Earth orientation parameters, GPS car- 
rier phase biases, and random-walk zenith troposphere de- 
lays for each tracking site; all transmitter/receiver clocks 
but one were treated as white noise parameters. X and 
Y polar motion, pole rate, and UT1-UTC rate are esti- 
mated as constant parameters (reset every 24 hr). 


The GPS constellation consists of satellites in Earth or- 
bit configured in six evenly spaced orbit planes. At times, 
satellites in certain orbit planes experience what will be re- 
ferred to in this article as eclipsing, or shadowing. Figure 2 
shows a satellite in an eclipsing orbit plane. GPS satel- 
lites eclipsing during GPS week 660 were PRN02, PRN14, 
PRN16, PRN20, PRN21, and PRN23. 

The results presented here illustrate recent improve- 
ments in GPS orbit accuracy. The three-dimensional 
root-sum-square ( SDRSS) orbit repeatability for each GPS 
satellite ( j ) is defined as [7] 


3DRSS(j) = 






<= 1 


( 1 ) 


One of the most important recent innovations is a new 
approach to modeling the effects of solar radiation pres- 
sure on the satellite orbits. The solar radiation environ- 
ment of the GPS satellites is generally constant except 
when a satellite’s orbit is in eclipse season. When this 
occurs, the satellites pass through Earth’s shadow, chang- 
ing the amount of solar radiation the satellite receives. In 
general, three body-fixed solar-radiation pressure param- 
eters are estimated for all GPS satellite orbits, regardless 
of whether or not those satellites are in eclipsing orbits. 
For this analysis however, the GPS orbits were estimated 
with five solar pressure parameters, which are shown in 
Table 2. Two solar radiation pressure parameters were es- 
timated as constant: GY C and GXjGZ , where X , Y, and 
Z represent spacecraft body-fixed coordinates and GY C 
is the “Y-bias” parameter [4]. GXjGZ represents a sin- 
gle combined scale factor for the ROCK4 2 solar radiation 
parameter [4], while GX and GZ are independent scale 
factors for the X and Z directions. The three remain- 
ing solar pressure parameters are estimated as stochas- 
tic corrections to the constant solar pressure parameters, 
modeled as a first-order Gauss-Markov process [5]. This 


where j is the satellite PRN number, T is the number 
of epochs within both overlap segments, t is the time in- 
dex, k indicates which overlap segment, and d? )ki is the 
three-dimensional distance between corresponding points 
on the two overlapping segments. In other words, the orbit 
quality of a single day is quantified as the rms difference 
between the ephemerides computed over the correspond- 
ing 3 hr of orbit overlap at both ends of that day. This 
concept is illustrated in Fig. 3. 3 

Figure 4 shows the SDRSS orbit repeatability for all 
GPS satellites active during GPS week 660. The improve- 
ment in the orbits with stochastic solar pressure param- 
eters is approximately 25 percent overall for all the GPS 
satellites (averaged), with a 44 percent improvement in the 
eclipsing orbits and only an 8 percent improvement in the 
noneclipsing orbits. This demonstrates how mismodeling 
of satellite force parameters due to solar radiation can have 
a significant effect on GPS orbit accuracy, especially for an 
eclipsing satellite. In an eclipsing orbit, a satellite does not 
experience constant solar radiation forces throughout the 
entire orbit arc. 


1 S. C. Wu, Y. Bar-Sever, S. Bassiri, W. I. Bertiger, G. A. Hajj, 
S. M. Lichten, R. P. Malla, B. K. Trinkle, and J. T. Wu, TOPEX/ 
POSEIDON Project: Global Positioning System (GPS) Precison 
Orbit Determination (POD) Software Design , JPL D-7275 (inter- 
nal document), Jet Propulsion Laboratory, Pasadena, California, 
March 1990. 

2 ROCK4 is a solar radiation pressure model provided by Rockwell 

Corp. for GPS satellites. 


The orbit repeatability for each satellite using both 
strategies is shown in Table 3. These values represent 


3 J. F. Zumberge, “Regular Analysis of Global Positioning System 
Data from a Globally- Distributed Network of Rogue Receivers,” 
Interoffice Memorandum 335.4-92-012 (internal document), Jet 
Propulsion Laboratory, Pasadena, California. 
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the average orbit repeatability over the 7 days in GPS 
week 660, with the eclipsing satellites indicated. This ta- 
ble shows the improvement in orbit accuracy due to esti- 
mating stochastic corrections to the GPS solar radiation 
parameters (srp). The table also shows how the eclipsing 
satellites benefit more from improved estimation strategy. 


IV. Concluding Remarks 

This analysis shows that GPS orbit accuracy can be 
improved by estimating stochastic corrections to the GPS 
dynamic parameters. This new approach improves orbit 
quality for eclipsing satellites from 85 to 47 cm. This level 
of orbit accuracy is in agreement with results given by 
Zumberge et al. [7], where routine processing of GPS data 
shows orbit accuracy in the range of 25 to 50 cm. A di- 
rect result of the improvements in orbit accuracy can be 
seen in the improvement of the geocentric station coor- 
dinate accuracy [6]. The goal of the analysis described 
in [6] was to achieve centimeter-level accuracy for global 
geocentric coordinates. Those GPS results were obtained 
with 3 months of GPS measurements and compared to 
satellite laser ranging (SLR) solutions from many years 
of repeated observations. It was demonstrated that the 
geocenter estimates from GPS are accurate to better than 
2 cm in the X and Y components and to approximately 


8 cm in the Z (where Z is parallel to the axis of rotation). 
This capability has important benefits for DSN tracking 
and for geophysical research, such as geocentric crustal 
motion studies, and for understanding the magnitude and 
time scale of geocenter variations and their origin. Precise 
tracking of interplanetary spacecraft and Earth orbiters 
requires that DSN geocentric station coordinates be de- 
termined to high accuracy. The accuracy for geocentric 
DSN coordinates determined from the GPS data is ap- 
proximately 3 to 10 nrad in angular measure. 

Figure 5 shows a history of the improvements in the 
GPS determination of the geocenter. Most of the re- 
cent improvements can be attributed to the changes in 
the technique used for modeling solar radiation pressure, 
changes that have been described in this article. These 
new results enable the tracking sites to be precisely spec- 
ified in a reference frame whose origin is at the geocen- 
ter, and will enable precise alignment of different reference 
frames used for Earth-based tracking, interplanetary navi- 
gation, and geophysical measurements. Some future work 
in this area will examine systematic error sources, further 
improve the orbit modeling, study fixing phase ambigui- 
ties, and incorporate data from low Earth orbiters, such as 
TOPEX/POSEIDON, to look for time-varying signatures 
of the geocenter. 
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Table 1. GPS Rogue receiver sites. 


ALBH 

Albert Head, B.C., Canada 

NALL 

Ny Alesund, Norway 

ALGO 

Algonquin, Ontario, Canada 

ONSA 

Onsala, Sweden 

CANB 

Canberra, Australia 

PAMA 

Pamatai, Tahiti 

FAIR 

Fairbanks, Alaska (USA) 

PENT 

Penticton, B.C., Canada 

GOLD 

Goldstone, California (USA) 

PGC1 

Victoria, Canada 

HART 

Hartebeesthoek, South Africa 

PINY 

Pinyon, California (USA) 

HARV 

Harvest Platform, California (USA) 

QUIN 

Quincy, California (USA) 

HERS 

Herstmonceux, Great Britain 

RCM2 

Richmond, Florida (USA) 

HONE 

Honefoss, Norway 

SANT 

Santiago, Chile 

JPLM 

Pasadena, California (USA) 

SCRI 

La Jolla, California (USA) 

KOKB 

Kokee, Hawaii (USA) 

ST JO 

St. Johns, Canada 

KOSG 

Kootwijk, The Netherlands 

TAIW 

Taiwan 

MADR 

Madrid, Spain 

TROM 

Tromso, Norway 

MASP 

Maspalomas, Grand Canary Is., Africa 

USUD 

Usuda, Japan 

MATE 

Matera, Italy 

WETB 

Wettzell, Germany 

MCMU 

McMurdo Station, Ross Is., Antarctica 

YAR1 

Yarragadee, Australia 

METS 

Metsahovi, Finland 

YELL 

Yellowknife, Canada 


Table 3. 3DRSS orbit repeatability for GPS week 660. 


PRN number 


With stochastic 
srp, m 


Without stochastic 
srp, m 


Table 2. Estimated solar radiation parameters. 


Parameter 

Model 

A priori sigma 

GX/GZ 

Constant 

100% 

GYC 

Constant 1 

2 nm/sec 2 

GX 

First-order G-M b 

10% 

GY 

First-order G-M b 

0.1 nm/sec 2 

GZ 

First-order G-M b 

10% 


a Estimated as constant parameter with no process 
noise. 

b Gauss-Markov (G-M) with time constant of 4 hr and 
steady-state sigma of 0.1 nm/sec 2 . Typical magni- 
tude of GX and GZ accelerations is 100 nm/sec 2 . 


2 a 

3 

11 

12 

13 

14 a 

15 

16 a 

17 

18 
19 
20 a 
21 a 
23 a 

24 

25 

26 
28 


0.76 

0.70 

0.38 

0.43 

0.32 

0.48 

0.51 

0.28 

0.52 

0.43 

0.66 

0.72 

0.30 

0.29 

0.45 

0.42 

0.30 

0.48 


0.64 

0.28 

0.32 

0.45 

0.56 

0.94 

0.58 

1.02 

0.56 

0.60 

0.44 

1.19 

0.75 

0.54 

0.56 

0.47 

0.83 

0.43 


a Eclipsing satellites. 
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Fig. 1. 1992 GPS tracking sites. 
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Fig. 2. Satellite in Earth orbit with simple cylindrical shadow model. 
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Fig. 3. Assessment of GPS orbit accuracy. 
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Fig. 4. 3DRSS GPS orbit repeatability— 7 days, GPS week 660. 



Fig. 5. Recent improvements in GPS geocentric coordinate 
accuracy. 
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